Selective Efficient Quantum Process Tomography 
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We present a new method for quantum process tomography. The method enables us to efficiently 
estimate, with fixed precision, any of the parameters characterizing a quantum channel. It is selective 
since one can choose to estimate the value of any specific set of matrix elements of the super-operator 
describing the channel. Also, we show how to efficiently estimate all the average survival probabilities 
associated with the channel (i.e., all the diagonal elements of its x~niatrix). 
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The efBcient characterization of the temporal evolution 
of a quantum system is one of the main tasks one needs to 
accomplish to achieve quantum information processing. 
In particular, this is essential to determine the most im- 
portant sources of errors to design appropriate quantum 
error correction techniques. The set of methods used to 
determine the evolution of a quantum system are generi- 
cally denoted as quantum process tomography (QPT) . In 
general, QPT is a very hard task whose completion typ- 
ically requires resources scaling exponentially with the 
number of qubits in the system (n). To show that this is 
indeed the case, we can proceed as follows: Under very 
general assumptions the evolution of the quantum state 
of a system can be represented by a linear, completely 
positive, trace preserving map pout — £{pin)- Choosing 
a base of operators {Em}, the map can be written as: 
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Thus, the map is completely characterized by the positive 
hermitian matrix Xmn that must also satisfy the above 
trace perserving condition. Therefore, fully character- 
izing the map requires — Z?^ real parameters (where 
D = 2" is the dimension of the Hilbert space of n qubits). 
For this reason QPT is a highly inefficient task. 

In this letter we present a method that enables us to 
evaluate (with fixed precision) any of the coefficients of 
the x-matrix with resources that scale polynomially with 
the number of qubits. The method is "selective" since 
one can use it to estimate any coefficient (or any set of 
coefficients). For each coefficient there is an efficient es- 
timation strategy that we describe below. For this rea- 
son we denote our strategy as selective efficient quantum 
process tomography (SEQPT). Our method is inspired 
on previous proposals that use randomized subroutines 
as intermediate steps for efficiently estimating any aver- 
age gate fidelity. Moreover, we also propose a method 
to estimate all the diagonal Xmm coefficients using only 
polynomial resources. 

First, it is convenient to review the main tomographic 
schemes and their properties. Standard quantum process 



tomography (SQPT) was the first tomographic method 
proposed [l|] . It involves preparing a set of input states 
Pk, and then performing full quantum state tomogra- 
phy on the resulting output states obtained after evo- 
lution. By doing this, we directly measure coefficients 
Xjk = Tr{pkS{pj)). However, if one wants to find the 
matrix elements Xmn it is necessary to invert an expo- 
nentially large system of equations relating A with x Wi ■ 
For this reason, the method is indirect (since it requires 
inversion to obtain matrix elements Xmn)- It is also in- 
efficient since, in the most general case, in order to esti- 
mate any of the coefficients Xmn one needs to perform an 
exponentially large number of experiments and classical 
postprocessing. 

Direct Characterization of Quantum Dynamics 
(DCQD) was recently proposed 0, |j| and it requires an 
extra ancillary system of n clean qubits with a clean 
quantum channel. Provided such a resource is available, 
the method enables the direct estimation of all diagonal 
Xmm by associating each result of a single complete 
measurement to each one. The method requires a rather 
small number of elementary quantum gates (0(n)). A 
prescription for estimating any Xmn is given. However, 
estimating off-diagonal coefficients may require the 
inversion of a system of equations which, for some 
off-diagonal Xmn coefficients, involves an exponential 
number of experiments. In these cases this method is 
inefficient. 

Symmetrized Characterization of Noisy Quantum Pro- 
cesses Q (SCNQP) transforms the channel E into a sym- 
metrized channel £' via twirling operations. After sym- 
metrization, only diagonal x'mm coefficients remain, be- 
ing the averages over the original coefficients of the same 
Hamming weight. The twirling is achieved using only 
(0(n)) single qubit gates with constant depth. The val- 
ues of the averaged coefficients are linearly related to 
output probabilities through an upper diagonal square 
matrix of size n -I- 1. The method is ideally tailored 
for evaluating the applicability of relevant quantum error 
correcting codes [i| . as it allows the evaluation of diago- 
nal Xmm coefficients averaged over operators of the same 
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Hamming weight (i.e. xooj average over 1 qubit errors, 
etc). However, it is not possible to estimate any of the 
off-diagonal Xmn coefficients, which are wiped out by the 
symmetrization protocol, nor distinguish among specific 
Pauli errors of the same Hamming weight. 

Our method has a similar flavor to SCNQP adding 
the possibility to determine any of the coefficients Xmn 
(including off-diagonal ones) with polynomial resources. 
The method is based on two observations: The first, is 
the fact that any matrix element Xmn can be related 
to an average survival probability of input states under 
the action of the channel (or a related quantity as de- 
scribed below). The average involved here is over the 
entire Hilbert space using the so-called Haar measure. 
The second observation is that such averages can be effi- 
ciently estimated by sampling over a finite set of states (a 
2-design, as described below). Let us describe these two 
crucial observations in more detail. Before doing that we 
point out that a choice of base to define the Xmn matrix is 
required. We will use an operator base {Em} satisfying 
X.r{EmEl) = D5m,n (orthonormality) and. E„,E}-^ = I 
(unitarity). For definiteness only, we assume Eq ~ I. A 
convenient choice is the base formed by tensor product 
of Pauli operators {Pm}, which have a convenient group 
structure: PmPn = i^™ "Pmitn, where ★ is the commuta- 
tive group operation. 

We start by noticing that the average fidelity of the 
map £ is defined as: [7| 
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In fact, F{£) is nothing but the survival probability av- 
eraged over all pure states \'4>)- ft is quite simple to see 
that there is a direct relation between the average fidelity 
F{£) and the coefficient xoo- As was shown in [8], for any 
operators Oi and O2 we have: 
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where P^ = lip) {ipl is the projector on state lip). Using 
this equation together with the x^matrix representation 
for the map £ in the definition of the average fidelity ([2]) 
we find F(£) — ^^l^^f^ ■ All other diagonal coefficients 
Xmm are directly related to average fidelities of slightly 
modified channels. In fact, if we consider the channel 
£m{p) = El^£{p)Em, its fidelity is given by: 



F{£m)= / {i^\El£{P^)E^\i;)di; 
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Note that all the above fidelities have a lower bound of 
1/{D + 1). To measure the diagonal coefficients Xmm we 
must perform the experiment described in Fig. f and 
average over all states |V')- 

The off-diagonal elements Xmn can also be related to 



FIG. 1: Circuit for measuring Xmm for a given channel £ 



average quantities by using the following identity: 



\£{ElP^E„)\^)d^Jj 
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This equation can be obtained by using the x~ 
representation for the map £, together with equation ([3|) 
and the trace preserving condition for £: 

J2 Xm'nMErn'ElEnEl) = tr{ElEn) = DSm,n (6) 
m'n' 

To measure the complex off-diagonal coefficients Xmn we 
can proceed as follows: We add an extra clean qubit and 
consider the map to be the one described by the cir- 
cuit in the dashed box of Fig. 2. Thus, if C-E^ denotes 
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FIG. 2: Circuit for measuring Re(xmri) for a given channel £ 

the controlled Em operation, we have £mn{(^) = £{C- 
E\C-El^HaHC-EmC-En) (where a is the joint state of 
the ancillary qubit and the original system). After ap- 
plying this operation we can use equation ([5]), and realize 
that the real part of Xmn can be obtained by measuring 
the polarization of the ancillary qubit (i.e. the expecta- 
tion value of (Ta; conditioned on the survival of the state 
■0}, averaging over all states Thus, this is due to 

the fact that 



tr{£mn (|0) (0| ® P^) a^ (g)P^)diP = 



-DRe(xmn) + S 

71 



D + l 

(7) 

If instead of measuring the expectation of (Jx P^p we 
measure the expectation value of ay ® P^, the average 
yields ^inix^^Hi) 



over 



D+l 



Summarizing: So far we showed that all matrix ele- 
ments Xmn can be directly related to average survival 
probabilities. In fact, the diagonal coefficients Xmm are 
related to average fidelities of the channel £m- The off 
diagonal coefficients are related to the average polariza- 
tion of the ancillary qubit conditioned to the survival of 
the state after the evolution with the channel £mn- 

Experimentally measuring these averages over the 
Hilbert space seems completely unrealistic. However, 
the beautiful recent work on the theory of 2-designs 
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^ 3 U , 3 provides the means for doing so. P. Delsarte 



13l | showed how integrating polynomials on the sphere 
could be reduced to averaging the integrand on a finite 
set of points coined spherical designs. The same idea can 
be extended to integrals over the entire Hilbert space. 
Here, we only need a state 2-design K, that satisfice 

/ (V-l OiFv02 IV-) d^^^\Y. (V'l OiP^02 IV) , (8) 

for all operators Oi,2- Thus, averaging over the entire 
Hilbert space is equivalent to averaging over the finite set 
X. State 2~design with a finite (but exponentially large) 
number of states exist. Although the computation of 
the exact average over the 2-design is still exponentially 
hard, it is now possible to realize that an estimate for 
the average can be efficiently found. This average can be 
estimated by randomly sampling over initial states 
chosen from the set X. This is the final piece of our 
method. 

Luckily, it is rather easy to produce a state 2-design 
for n qubits. One possibility is finding D -\- 1 mutually 
unbiased bases (MUB) that automatically form a state 
2-design Each base will be labeled with an index 

J = 0, . . . , and the states within each base will be 
labeled with the index m = \, . . . , D. In order for the or- 
thonormal bases to be unbiased, the D{D + 1) states of 
the MUBs must satisfy = ^ for all J K. 

Since generalized Pauli operators may be partitioned into 
D + 1 maximally sets of D commuting operators so that 
each pair of sets only hold the identity / as common ele- 
ment there are D + 1 MUBs, each one diagonalizing 
each of these commuting subsets of Pauli operators ■ 
In this way, the set of states in the MUBs can be ef- 
ficiently described and also can be efficiently generated 
with O(n^) one and two qubit gates It is simple 

to adapt the procedure used to efficiently generate any 
state in any MUB to compute survival probabilities of 
such states and also to compute the transition rates from 
the (J, m) to (J, m') states. 

Other strategies may use 2-designs other than MUBs. 
For example, Dankert et. al. [1, IlBl propose to use of 
approximate unitary 2-designs (which are designs on the 
space of unitary operators) showing that they can be effi- 
ciently approximated. An approximate unitary 2-design 
with accuracy e -I- l/D^ can be obtained by employing 
only 0{n log i) gates. Unitary 2-designs acting on any 
fixed state induce state 2-designs fitting into the previous 
scheme. Dually the action of the random unitaries may 
be interpreted as symmetrizing the channel £ through 
twirling. Following this line, we may also use weaker 
symmetrization protocols as in SCNQP for estimat- 
ing fidelities of modified channels (|ll2p . 

Let us summarize the complete method. The estima- 
tion of diagonal coefficients Xmm of the quantum map £ 
can be efficiently done by evaluating the fidelity of the 



channel £m averaged over a random sample of the 2- 
design formed by the states of D -I- 1 MUBs. The real 
and imaginary parts of the off-diagonal elements Xmn are 
obtained from the x and y-polarizations of the ancillary 
qubit conditioned on the survival of the states randomly 
chosen from the same 2-design when the evolution of the 
combined system-|-ancilla is described by the map £mn- 
The strongest requirements for this method to be applica- 
ble are 2-design state generation and the implementation 
of the single qubit controlled Pauli operators used in the 
construction of the £mn channels. 

The fact that the estimate of Xmn can be obtained by 
sampling a small subset of states is crucial for the effi- 
ciency of the method. For definiteness, suppose we are 
interested in estimating -^''^(xrnre)+i ^^-^^Yi precision e. The 
results of a given experimental run may produce three 
outputs: 0, corresponding to no survival on the main 
system, and ±1 corresponding to survival and different 
polarizations eigenstates on the auxiliary output. Hence, 
independent measures have a variance smaller than 1. 
Averaging over M > independent experiments, we 
may obtain the desired precision. Similarly the estima- 
tion of fidelities F{£m) may achieve a precision e using 
only M > independent experiments. The impor- 

tant result is that the number of required experiments 
does not increase with D. 

Up to this point we have only considered measuring 
survival probabilities. However, after preparing a state 
l^fc)' '■^^^ measure transition probabilities to other 
states I'f/'fc')- We will show that in this way we can effi- 
ciently estimate all the diagonal coefficients Xmm- To do 
this, it is crucial to use the 2-design formed by the MUBs 
associated to the same operators used in the x~rnatrix 
description of £ (the Pauli operators are well suited for 
that purpose). Suppose that we perform M experiments 
and we store the results as triplets (J, k, k') correspond- 
ing to the labels of the initial state and the one detected. 
We will show that the fidelity F{£jn) can be estimated 
as the frequency of certain events (events of the m-type). 
Below, we will describe a simple procedure to decide if 
a given event is of the m-type or not (the procedure re- 
quires resources polynomial in n). 

Each base is determined by n commuting Pauli opera- 
tors Pi . . . which are the generators of the Abelian 
group of operators that are diagonal in that particu- 
lar basis. We can label the base elements {\ip'l)) ac- 
cording to the eigenvalues of the n generators. Thus, 
the label A: is an n components binary vector: the i- 
th component ki determines the ±1 eigenvalue of Pf as 
P/ = (— l)*^' IV'fc )• The action of any Pauh P on the 
state \tp^) transforms it into another state of the same 
base J. Those transition rules are fully determined by 
the binary vector p that encodes the commutation pat- 
tern between P and the generators of the J-th base (i.e. 
the vector p is such that (-I)P'PP/ = P/ P ). Thus, 
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'^k+p/ (^'^P ^ phase), where k and p are 
added bitwise modulo 2. Hence, the probability to detect 



the state 



i+p) 



in the final measurement after evolving 



with any unitary operator U is identical to the proba- 
bility of detecting the state 1-0^) after evolving with the 
operator obtained as the product P x U. Applying this 
observation we conclude that in order to measure any 
diagonal coefficient Xmm we could modify the strategy 
we described above by removing the final Em gate pro- 
ceeding as follows: We compute Xmm by averaging the 

whenever the 



k+Pr, 



transition probability to the state 

state is prepared {pm is the commutation vector of 
Em and the generators of the base J) . 

The procedure to estimate all diagonal Xmm coeffi- 
cients must be clear now. To estimate the associated 
fidelity F{£m) from the experimental data obtained af- 
ter M experiments we must compute the frequency of 
the events (J, k, k') that satisfy the condition k + k' = pm 
(where Pm is the commutation vector of the operator Em 
and the generators of the base J). These are the m-type 
events. Performing this check requires knowing the vec- 
tor Pm (which requires O(n^) classical operations). The 
lassical complexity for this process is 0{n'^M). However, 
to the best of our knowledge, the construction of the n 
generators for base J requires 0{n'^) classical operations, 
which is the dominating overhead. Errors corresponding 
to the estimations of the different Xmm are correlated. 
However the variance of any such estimator behaves in 
the same way as the one corresponding to the method to 
evaluate a single Xmm at a time. 

We can also use this method to devise a test to effi- 
ciently detecting the coefficients Xmm with values that 
are above a certain threshold (which is D -independent, 
i.e. not exponentially small). The idea is that if the co- 
efficient Xmm Is large then the condition k + k' = Pm 
must be satisfied frequently in the experimental data. 
Given the triplets corresponding to two experiments per- 
formed on different bases, we may efficiently determine 
the unique operator Em satisfying the condition for both 
of them (using O(n^) operations). The number of such 
pairs of triplets is bounded by M{M + l)/2. Therefore 
we can find out the operators Em satisfying the condi- 
tion k + k' = Pm for at least two bases J with an over- 
head that depends at most quadratically M. All the 
operators satisfying such criterion will be obtained with 
0{M'^n^) classical operations. After sieving the opera- 
tors in this way we can focus on those that passed the 
test to study them further. The conclusion is that in this 
way, it is possible to efficiently identify and estimate all 
large Xmm coefficients. Again, we stress that this method 
enables us to estimate if the coefficients are larger than a 
fixed D-independent threshold such as M~^/^. We may 
also simultaneously estimate groups of off diagonal Xmn 
coefficients for which E^En is the same. The required 
book-keeping is similar to the one used for diagonal co- 



efficients but requires additional care regarding the nu- 
merical phases that arise. 

Outlook- We have shown how any particular Xmn coef- 
ficient for a channel £ may be estimated from the average 
survival probabilities of at most two modified channels 
requiring at most one extra clean qubit. This approach, 
together with the use of randomized fidelity measurement 
schemes is the first to allow the efficient estimation of any 
desired Xmn coefficient with a number of experiments 
that is dimension independent. The method shows its 
strength when only partial tomography is desired, since 
the resources required are polynomial in the number of 
subsystems. We further presented a possible extension by 
sampling over initial states belonging to the MUBs asso- 
ciated with the same operator base used in the channel 
representation. This allows us to profit from the informa- 
tion provided by complete base measurements. Finally, 
we prove the possibility of efficiently detecting and char- 
acterizing Pauli channels with sparse x-matrices. JPP 
is a Fellow of CONICET. This work was supported with 
grants from Anpcyt and Conicet. 
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